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Abstract 

We report predominantly zero field muon spin relaxation measurements in 
a series of Ca-doped LaMnC>3 compounds which includes the colossal mag- 
netoresistive manganites. Our principal result is a systematic study of the 
spin-lattice relaxation rates 1/Ti and magnetic order parameters in the se- 
ries Lai-zCazMnC^x = 0.0,0.06,0.18,0.33,0.67 and 1.0. In LaMn0 3 and 
CaMnOs we find very narrow critical regions near the Neel temperatures 
T/v and temperature independent 1/Ti values above Tat. From the 1/Ti in 
LaMnC>3 we derive an exchange integral J = 0.83 meV which is consistent 
with the mean field expression for T\r. All of the doped manganites ex- 
cept CaMnOs display anomalously slow, spatially inhomogeneous spin-lattice 
relaxation below their ordering temperatures. In the ferromagnetic (FM) in- 
sulating Lao.82Cao.i8Mn03 and ferromagnetic conducting Lao.67Cao.33Mn03 
systems we show that there exists a bi-modal distribution of /xSR rates A/ 
and A s associated with relatively 'fast' and 'slow' Mn fluctuation rates, re- 
spectively. A physical picture is hypothesized for these FM phases in which 
the fast Mn rates are due to overdamped spin waves characteristic of a disor- 
dered FM, and the slower Mn relaxation rates derive from distinct, relatively 
insulating regions in the sample. Finally, likely muon sites are identified, and 
evidence for muon diffusion in these materials is discussed. 

PACS numbers: 71.27.+a, 72.15.Qm, 76.75.+i, 75.30.Mb 
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I. INTRODUCTION 

The strong coupling of the spin, charge and lattice degrees of freedom in the colossal 
magnetoresistive (CMR) manganites gives rise to a wide variation in ground-state behavior, 
including paramagnetism, ferro-and antiferromagnetism, charge ordering (CQ1, and insulat- 
ing and conducting charge transport]!!. The earliest studies of these systems™ involved the 
substitution of a divalent alkaline earth (R = Sr, Ca, Ba) for trivalent La in the perovskite 
structured LaMnOs, an insulating antiferromagnet. This substitution produces a hole at a 
Mn site (i.e., Mn 3+ changes to Mn 4+ ) which can hop between Mn ions with different valence 
states via the ferromagnetic double exchange (DE) interactional. Linkage between the DE in- 
teration and charge mobility provided an early explanation of the paramagnetic/insulating to 
ferromagnetic/conducting transition in Lai_ x R x Mn03. More recent experiments^ and sub- 
sequent theoretical analyses^ have also established the importance of the local Jahn- Teller 
(JT) structural perturbation in the Mn 3+ state. Thus, both the JT and DE interactions are 
important in determining the ground state properties and magnetic ordering temperatures 
in these materials. The JT and DE interactions are also thought to induce the formation of 
magnetoelastic polarons, local structural distortions surrounded by polarized Mn spins, in 
CMR compounds. The degree of polaron formation depends upon the degree of cation size 
mismatch between the La atom and the alkaline earth atom.Q Sr produces the least distor- 
tion and Ba the most. These concepts comprise the fundamental ideas behind the strong 
coupling of spin, charge and lattice degrees of freedom in CMR materials. In this paper 
we study the spin dynamics of the Lai_ x Ca x Mn03 system using the muon spin relaxation 
(//SR) techniqueja thus providing a valuable complement to other experiments probing the 
charge and lattice degrees of freedom in the these materials. 

The Lai_ z Ca x Mn03 system has been extensively studied because it exhibits a wide 
variety of ground states as a function of x below room temperature B At temperatures be- 
low about 700 K the system is orthorhombic for all Ca concentrations. The end members 
LaMnOs and CaMn03 are both antiferromagnetic (AFM) insulators with different mag- 
netic structures. The addition of Ca to LaMnOs produces a canted AFM insulating state 
for < x ;$ 0.07, a ferromagnetic (FM) insulating state (followed by CO at lower tem- 
peratures) for 0.07 < x ;$ 0.21, and a FM metallic state for 0.21 < x ^ 0.50. Above 
x = 0.50 the system remains insulating with a high temperature CO transition observed 
for 0.50 < x < 0.875, followed by AFM at lower temperatures. Experimentally, transport 
measurement sEj are consistent with the formation of magnetoelastic polarons in the para- 
magnetic state of Lai_ x Ca x Mn03, and both neutron scattering^ and x-ray fine structure 
analysis^ find evidence for the persistence of JT distortions significantly below the FM tran- 
sition. Neutron scattering has been used to study the low-temperature magnetic properties 
in FM perovskites,E3 and has revealed a broad peak centered around zero energy transfer 
which coexists with Mn spin waves near TcJiil 



A. /tSR experiments and data analysis 

The /iSR technique employed here involves the intersititial implantation of positive 
muons with essentially 100% spin polarization, oriented antiparallel to the initial muon 
momentum vector. The /i + decays with a mean lifetime r^ of 2.2 /xs into a positron and two 



undetected neutrinos. The positron is emitted preferentially along the muon spin direction; 
detection of the positron emission rate in counters aligned along the initial muon spin di- 
rection allows one to monitor the time rate of decay of the muon spin polarization along its 
initial spin direction. The positron emission rate dN/dt is given by 

dN/dt = B + iV (l/r M ) exp(-W)[l ± AG z (t)], i 1 ) 

where B is a background (usually independent of time), N is a normalization factor, A is the 
average asymmetry of the decay angular distribution (typically 0.2—0.3), and G z (t) describes 
the time rate of decay of the muon spin polarization. The + or — signs are appropriate for 
positron counters in the direction of, or opposite to, the muon spin, respectively. 

The /zSR data presented here were taken at the M20 muon channel at TRIUMF (Van- 
couver, Canada) and on the GPS spectrometer at the Paul Scherrer Institute (Villigen, 
Switzerland). These data could be fit to a relaxation function G z {t) = G osc (t) +G T \ x (t), cor- 
responding to oscillating and relaxing terms, respectively. In zero applied field the oscillating 
component G osc (t) occurs in a magnetically ordered state and is given by 

G osc (t) = Sj(A osc )j exp(-t/T 2i ) cos(27rz/^t + 0^), (2) 



in 



where (A osc )i is the amplitude of the i precessing component, (j)^ is a phase angle and z/ ( 
and 1/T 2 j are the corresponding muon precession frequency and inhomogeneous damping 
rate. We define Sj(A osc )j = v4 osc . The existence of multiple components can be due to more 
that one muon lattice site and/or different microscopic magnetic environments (even for a 
single muon lattice site). The quantity litv^ = j^Bi, with 7 M the muon gyromagnetic ratio 
(8.54 x 10 4 Hz/Oe) and B t the local magnetic field. 

In our experiments the relaxing component G r \^(t) frequently approximates a stretched- 
exponential form: 

G r ix(t) = A rlx exp[-(t/T 1 )^], (3) 

where 1/Ti is a characteristic spin-lattice relaxation rate. The polycrystalline averages of 
A r i x and A osc in a magnetically ordered state (where both A r \ x and A osc are non-zero) are 
A r \ x = 1/3 and A osc = 2/3. For rapid fluctuations the muon relaxation rate 1/Ti is given 
by 1/Ti oc 7^D(j \fiB(q)\ 2 T(q), where \5B(q)\ is the amplitude of the fluctuating local field 
and r{q) is the Mn-ion correlation time. The sum over the characteristic momenta q of 
the Mn-ion excitations occurs because the muon is a local probe. When the exponent 
K < 1 in Eq. (3) a distribution of 1/Ti values is implied,^ and thus |&B(g)| and/or r(q) 
are distributed. Typically, near a magnetic phase transition the correlation time becomes 
longer, causing 1/Ti to increase, a phenomenon known as critical slowing down. Below the 
ordering temperature, the fluctuating amplitude decreases and 1/Ti is reduced. 

The stretched exponential form in Eq. ([J) is often useful to parameterize the data 
because it involves only a few parameters and can be used to approximate G r ix(£) over a 
wide temperature range. However, for the materials under investigation here, we have found 
that for large 1/Ti and small K the relaxation function can often be better approximated 
over a limited temperature range by the sum of two exponentials: 

G t]x (t) = A f exp(-A / t) + A s exp(-A s t), (4) 



where A/ and X s correspond to fast and slow local-field fluctuation rates r _1 , respectively. 
This fitting function is for a bi-modal distribution of fluctuation rates, and its appropriate- 
ness is discussed in detail in Section III below. 



B. Sample preparation 

The Lai-^Ca^MnOs samples studied here were polycrystalline materials prepared ei- 
ther as boules in an optical floating zone furnace (0 < x < 0.18) or as sintered powders 
(0.33 < x < 1.0). In general, polycrystalline materials of (La,Ca)Mn03 are more composi- 
tionally homogeneous than comparable volumes of single crystals, because of relatively less 
Ca evaporation and more stable growth conditions. The samples were characterized using 
x-ray diffraction, resistivity and susceptibility. All samples were greater than 98% single 
phase. 



II. EXPERIMENTAL RESULTS: AN OVERVIEW 
A. LaMnOs and CaMnOs: Antiferromagnets 

The compounds LaMnOs and CaMn03 are the AFM end members of the (La,Ca)Mn03 
series. They are both insulating, but do not exhibit charge ordering. LaMnOs possesses 
A-type AFM ordera below the Neel temperature T/v ~ 139 K, in which ferromagnetic (FM) 
MnC>2 sheets are antiferromagnetically aligned with one another. Conversely, CaMnOs pos- 
sesses G-type AFM orderB below T/v ~ 123 K, where each Mn atom is antiferromagnetically 
aligned with its nearest neighbor. 

Figure |1] shows the temperature dependence of the two muon frequencies observed in 
LaMnOs in zero applied field, together with their amplitudes and fractional linewidths 
l/(27rz/ Ml T2j). These two frequencies have the same temperature dependence, and their 
signal amplitudes (A osc )j are essentially the same and independent of temperature. Note 
that the fractional linewidths below T/v are very narrow (< 2 %), as seen in the bottom 
frame of Fig. [I]. 

Figure [^ shows the temperature dependence of the spin-lattice relaxation rate 1/Ti and 
the relaxing amplitude A r \ x in LaMnOs. These data were obtained by fitting the relaxing part 
of G z (t) to Eq. (3). The exponent K was found to be consistent with a simple exponential 
relaxation function (K = 1) over the entire temperature region measured. The sharp rise 
in 1/Ti at T/v and fall just below T/v indicates a very narrow range of critical temperatures 
for this compound. The magnitude of 1/Ti is temperature independent above T/v. Note 
also that A r i x is equal to one above T/v, but falls within a few degrees to its polycrystalline 
average of about 1/3 just below T N , indicating that the transition to the AFM state is quite 
sharp in temperature. Both the rapid fall in A r i x and the sharp rise in 1/Ti occur very near 
the temperature where the observable muon frequencies approach zero. The latter is shown 
in the bottom frame of Fig. |2], where the two muon frequences have been normalized at 120 
K to show that they have the same temperature dependence. Finally, the top frame of Fig. |2| 
shows that a rather modest field of 3 kOe applied parallel to the muon spin polarization 
destroys the sharp critical behavior in 1/Tl. 



Figures ^ and ^ show comparable data for CaMnOs. In this compound three muon 
frequencies are observed at low temperatures, but only the lowest frequency line v^ is 
clearly observable at temperatures above about 90 K. The magnitude of v^ approaches zero 
frequency near T/v, as seen in the middle frame of Fig. |3| and the bottom frame of Fig. |j. 
The relative linewidths Xji^/KV^JT^i) are somewhat broader in CaMnOs than in LaMnOs. 
The dynamical signal in CaMnOs is shown in the top two frames of Fig. |], and exhibits 
a behavior similar to that in LaMnOs; e.g., a very narrow critical region, as indicated by 
sharp transitions in 1/Ti and v4 r i x , and a temperature-independent 1/Ti behavior above T/y. 

Summarizing our results on the two end compounds, LaMnOs and CaMnOs exhibit sharp 
transitions into their magnetically ordered states, and their dynamic relaxation functions 
G T i x (t) can be described by a single exponentially decaying time dependence. 



B. Lao.94Cao.o6Mn03: Insulating Canted Antiferromagnet 

As Ca is substituted for La one introduces chemical disorder into the lattice and holes 
into the Mn 3+ electronic structure, producing Mn 4+ . An additional consequence is that 
the local JT distortions are reduced as the holes become more abundant, because Mn 4+ is 
not a JT-active ion. In Lao.94Cao.o6MnOs, for example, the AFM order found in LaMnOs 
begins to show traces of FM, although only through a canting of the Mn spin momentsB 
The system remains insulating, although its resistivity drops in magnitude compared to the 
undoped LaMn0 3 (Fig. ||). 

The /xSR data for this system are now discussed. First, no oscillating component was 
observed. This is likely due to a complicated magnetic structure which washes out observa- 
tion of a discreet muon frequency or frequencies. The /iSR data for the relaxing component 
G r i x in Lao.94Cao.o6Mn03 are shown in Fig. |6|. These data were obtained by fitting G r i x with 
the stretched-exponential form in Eq. (||). One sees that the amplitude A r \ x drops sharply 
and the spin-lattice relaxation rate 1/Ti shows a peak at T « 122 K, indicating the onset of 
magnetic order. Note that Ca doping has decreased the transition temperature compared 
to the undoped LaMnOs. Comparing the zero-field relaxation rates in Lao.94Cao.o6MnOs to 
those in LaMnOs one notes that the width of the peak in 1/Ti is larger in the Ca-doped 
compound. This is so despite the fact that the transition width, as measured by the reduc- 
tion in the relaxing amplitude A v \x near the ordering temperature, remains quite narrow. 
Furthermore, the relaxation function in Lao.94Ca .o6MnOs is no longer exponential, as it was 
in LaMn0 3 . This is shown in the middle frame of Fig. ^, where the temperature depen- 
dence of the exponent K in Eq. (3) is displayed. G r ix actually begins to deviate from its 
exponential form (K = 1) at a temperature slightly above the peak temperature in 1/Ti. 
The value of K continues to decline as the temperature is lowered below the transition. As 
shown here and below, this non-exponential behavior and broadening of the transition in 
1/Ti is characteristic of the systems containing both La and Ca. Finally, as in both the 
undoped LaMn0 3 and CaMn0 3 compounds, the application of a small applied field tends 
to destroy the signature of the magnetic ordering in the spin-lattice relaxation rate (Fig. ||, 
bottom frame). 



C. Lao.82Cao.i8Mn03: Ferromagnetic Insulator 

The data in Lao^Cao.isMnOs, in which the level of doping is sufficient to induce ferro- 
magnetism, is now discussed. As shown in Fig. ^ the magnitude of the resistivity in this 
compound is lower than that of the Lao.94Cao.o6Mn03 material, indicating that the addition 
of doped holes is pushing the system closer to a conducting state. The material remains 
insulating, however. According to the phase diagram of Cheong and HwangB this material is 
a FM insulator with a transition to a CO state below about 60 K. Note that the resistivity 
for this compound shows a marked upturn at temperatures below about 80 K (see Fig. (|5])). 

Figure [7| shows the temperature dependence of v4 rlx , K and 1/Ti resulting from a fit to 
a stretched-exponential function for G v \ x . A sharp peak in the spin-lattice relaxation rate is 
evident at temperatures near 182 K, and a broader peak is evident at lower temperatures 
near 110 K. The sharpness of the peak at 182 K is characteristic of critical slowing down near 
a magnetic phase transition. This transition to an ordered state is also reflected in the sharp 
reduction in the magnitude of A r \ x below the same temperature. The peak temperature is 
consistent with the ferromagnetic transition temperature in the published phase diagram,0 
indicating that Tc = 182 K for this material. As in the Lao.94Cao.o6Mn03 compound, the 
exponent K decreases at temperatures below Tc- 

Figure |8| shows the temperature dependence of the measured muon frequency z/^, together 
with its fractional linewidth 1/ (2:7^^2 ) from the oscillating portion of G r i x . For comparison, 
the temperature dependence of 1/Ti, shown in the bottom panel of the same figure, exhibits a 
peak where v^ tends to zero. Note that only a single frequency is observed, unlike in LaMn03 
or CaMnOs, which display multiple muon frequencies. This indicates that the Ca doping is 
sufficient to have created a rather simple (FM) magnetic structure, so that a discrete muon 
frequency is observable, unlike the case of Lao.94Cao.o6Mn03. This observable frequency 
clearly corresponds to the buildup of the FM magnetization below Tc, but observation of 
the FM order parameter is only visible in the temperature range 150 < T < 180 K. At 
the lowest temperatures where v^ can be observed, T « 0.82Tc, the fractional spread in 
frequencies (l/(27rt / /i T 2 )) is at least an order of magnitude larger than found in the undoped 
LaMnOs. This presumably reflects the local disorder in the material, which shows up as a 
spread in local magnetic fields, even though the transition temperature itself remains fairly 
sharp. 

Below 150 K the ability to resolve a discrete muon frequency in Lao^Cao.isMnOs van- 
ishes, possibly due to an increasing linewidth or to a change in the magnetic state. At still 
lower temperatures, the muon relaxation rate 1/Ti shows a broad maximum, indicating a 
change in the Mn-ion spin dynamics. To observe this broad maximum it was necessary to 
freeze the value of K in the least-squares fitting procedure so that the form of the relaxation 
function would remain the same over the temperature range below about 170 K. 

D. Lao.67Cao.33Mn03: Ferromagnetic Conductor 

For Lao.67Cao.33Mn03 the Ca concentration is sufficient to produce a conducting state 
below the FM transition temperature Tc ~ 270 K (see Fig. |5|). These data are from the same 
sample as used previously;li3 however, the current data are of greater statistical precision 
than those previously reported and were taken at smaller temperature intervals near Tc- 



Figures |9] and |10] show the parameters obtained from fitting the relaxation function G z (t) 
using the stretched exponential form in Eq. (3) for G r i x - Only a single muon frequency v^ 
is found, the temperature dependence of which is shown in the middle frame of Fig. |9|. The 
fractional inhomogeneous linewidth 1/(270^2) is independent of temperature and yields 
about the same fractional width as in the insulating FM compound Lao^Cao.isMnOs. The 



spin-lattice relaxation rate 1/Ti is plotted in Fig. |10|, together with the exponent K and the 
amplitude A v \ x . The latter shows that the FM transition is not as sharp in this material as in 
the previously discussed compounds. Furthermore, the exponent K falls dramatically below 
T c , reaching values as low as about 0.2. Finally, the temperature dependence of 1/Ti reveals 
a rather broad peak, with large relaxation rates extending over a temperature range of nearly 
75 K below Tc- This persistence of a sizeable relaxation rate below the magnetic ordering 
temperature is more prominent at this Ca concentration than at any other presented here. 
Finally, a small applied magnetic field of only 1 — 3 kOe destroys the peak in 1/Ti near Tc- 

E. Lao.33Cao.67Mn03: Charge Ordering and Antiferromagnetism 

When 0.50 < x < 0.88 in Lai-^Ca^MnOa the system undergoes a CO transition, fol- 
lowed by AF order at lower temperaturesB The highest CO transition temperature appears 
to occur for x ~ 2/3 where T co ss 270 K and T/v ~ 125 — 150 K. These numbers are de- 
termined largely from dc susceptibility and resistivity measurement sB The /iSR relaxation 
data between 5 < T < 250 K for Lao.33Cao.67Mn03 are shown in Fig. [11], again using a 
stretched exponential form for G r ix- No /iSR oscillations were observed in this material. 
The temperature dependence of the spin-lattice relaxation rate 1/Ti shows two peaks, one 
near the Neel temperature T/v ~ 150 K and a lower, broader peak centered around 80 K. 
The temperature dependence of the relaxing asymmetry A T \ X shows a decrease from 100% to 
about 33% near 150K, which is the signature of magnetic order in a polycrystalline material. 
This onset of magnetic order coincides with the peak in 1/Ti, but the decrease in A rlx occurs 
over a 20 — 30 K temperature range, which is comparable to that in Lao.67Cao.33Mn03, but 
considerably larger than in the lower-doped samples discussed above. Also, similar to the 
observations in the previously discussed Ca-doped compounds, the exponent K decreases 
from about 1 above T/v to K < 0.4 below T/v. 



F. Ferromagnetic Materials Lao.82Cao.i8Mn03 and Lao.67Cao.33Mn03: A 

two- exponential analysis 

As mentioned above, the stretched-exponential relaxation function is characterized by 
only a few parameters and can therefore provide a useful parameterization of the data over 
a broad temperature range. For example, a deviation from the expectedE-3 and usually 
observedlia simple exonential form (K = 1) is readily detectable. However, careful exami- 
nation of the data in Figs. [| and [10] raises some concerns about the particular applicability 
of this function in Lao.67Cao.33Mn03. First, the exponent K reaches a very small value just 
below Tc, where K < 0.2. Even in spin glasses such as AgMn, where a stretched expo- 
nential /iSR relaxation function is observed, the exponent 1/3 < K < 1x3 (In AgMn a 
distribution of energy barriers for Mn spin rotation yields a broad distribution of correlation 



times r.) Second, in Lao.67Cao.33Mn03 the peak in 1/Ti actually occurs at a temperature 
below Tq (Figs. ||] and 1C), rather than being at Tc as expected. This effect is not seen in 
the other materials and most likely occurs because, in the temperature region where K is 
rapidly changing, the 1/Ti values at different temperatures are not derived from the same 
functional form. These issues suggest that a more refined analysis may be useful. 

Under ideal conditions determination of the precise shape of a decaying curve requires 
observation of the decay function over many decades with very high statistics.^ Although 
/iSR data do not always meet these stringent criteria, we find that the relaxation function G r i x 
can be better approximated by a sum of two exponentials (Eqn. 4), rather than a stretched 
exponential, when K < 1 and 1/Ti is reasonably large (usually > 0.1/xs -1 ). These conditions 
are realized for the two ferromagnetic materials Lao^Cao.isMnOa and Lao.67Cao.33Mn03, 
over a limited temperature range. This is illustrated in Fig. [12|, which displays G T \^(t) in 
Lao.67Cao.33Mn03. The top frame shows least squares fits for a simple exponential (K = 1) 
and a stretched exponential (K = 0.63). Clearly the stretched exponential function yields 
a better fit to the data. The bottom frame in this figure shows a fit to the sum of two 
exponentials using Eqn. 4; one can see at early times that the two-exponential fit is slightly 
better than that using the stretched exponential form. 

To more precisely illustrate the two-exponential nature of G r ix in La .67Cao.33Mn0 3 near 
Tc, in Fig. [13] we plot — t/ln(G r i x ) versus t on a log-log plot, together with the calculated 
curves from the fits to the stretched- and two-exponential functions shown in Fig. ^. Plotted 
in this way, a purely exponential relaxation function yields a horizontal line (zero slope); 
as the exponent K of the stretched-exponential is decreased from one, the calculated curve 
remains linear with slope (1— K). Clearly, the two-exponential model function yields a better 
representation of the data. An analysis of the relaxation function in Lao^Cao.isMnOs yields 
similar results. 

As stated above, fitting a relaxing curve to the sum of two or more exponentials is not 
rigorously unique. Credible results are obtained for Ca concentrations x = 0.18 and 0.33 
when K < 1 and 1/Ti > 0.1/is _1 , however, because the two rates A/ and X s in Eq. 4 



differ by more than a factor of ten. This is illustrated in Fig. [14]. For Lao^Cao.osMnOs 
and Lao.33Cao.67Mn03, although a stretched-exponential relaxation function fits the data 
reasonably well, the rates 1/Ti were too small to yield good, convergent two-exponential fits 
over a temperature range wide enough to produce interpretable results. 



Figs. ^ and [L6| show the temperature dependence of the parameters A/ and X s obtained 
from fitting G r i x in La .82Cao.i8Mn0 3 and La 067 Cao.33Mn0 3 , respectively, using Eqs. 3 and 
4. The two-exponential fits could only be performed over a limited temperature range above 
and below Tc, where the relaxation rate 1/Ti is large enough, as explained above. Sufficiently 
far above Tc, where K « 1, the data are well represented by a single exponential function. 
The open symbols correspond to the temperature region where two-exponential fits were 
feasible. The closed symbols at higher temperatures (top frames Figures [L^ and |16|) are 
from the single exponential fits. In order that all of the plotted data for Lao.67Cao.33Mn03 
in the top frame of Fig. [Tjd| correspond to an exponential functional form, some data near 
180 K were omitted from this plot, because a two-exponential fit was not feasible and the 
exponent K in Eq. 3 was significantly less than 1. The temperature dependence of the 
parameters Af and A s are shown in Fig. [17|. Here Af + A s = A r i x . 

Analyzed in this way the data for Lao^Cao.isMnOs and Lao.67Cao.33Mn03 yield very 



interesting results. First we note that in both materials only A/ shows a peak at Tc, 
characteristic of the slowing down of spins near their magnetic ordering temperature. The 
temperature dependence of X s shows no such peak, displaying only temperature indepen- 
dent behavior below Tc- Furthermore, the relaxation rates Xf and X s differ by at least an 
order of magnitude, with X s ~ 40A/ below Tc- The ratio of amplitudes Aj/A s near and 
below Tc is ~ 3.0 in Lao^Cao.isMnOs and ~ 1.5 in Lao.67Cao.33Mn03. Finally, as seen 
in Fig. [T7|, the amplitudes remain approximately independent of temperature below Tc in 
Lao.82Ca .i8Mn0 3 , whereas Af increases and A s decreases below T c in La .67Cao.33Mn0 3 . 

III. ANALYSIS AND DISCUSSION OF RESULTS 

A. Muon lattice sites 

The observation of two muon frequencies below T/v in LaMnOs indicates that the muon 
occupies magneticallv inequivalent sites in this material. A similar situation was observed 
by Holzschuh et al.fQ who used /iSR to study the muon's behavior in a series of single- 
crystal orthoferrite REFe03 (RE = rare earth) materials having the same crystal symmetry 
as the (La,A)Mn03 samples used in our experiments. Holzschuh et al. found that the muon 
occupies three different crystalline sites at low temperatures in both ErFeOs and YFe0 3 
but that, due to thermal activation, only one of these sites remains stable at temperatures 
above 100 and 300 K, respectively. This most predominant and stable site is located about 
1 A from the 0(1) oxygen atoms (as defined below) in the REO2 plane. 

As stated earlier the observed muon frequency v^ is given by (7^/271") .B, where the local 
magnetic field B arises from the sum of a transferred hyperfine field Bh yp and a local dipolar 
field B^p. In oxides, where the muon is bound about 1 A from an oxygen atom, the 
hyperfine field is in part due to the formation of a muon-oxygen covalent bond.El Holzschuh 
et al. found that their data in the orthoferrites were consistent with Bh yp <SL -£?<%, however. 
A similar finding has been reported from /iSR experiments in the cuprates.^ It is therefore 
reasonable to assume that Bh yp is also negligible in the manganites. 

Using the recently measured^ orthorhombic structure of LaMn03, where the Mn sit at 
(0,0.5,0) and the oxygen 0(1) at (0.073,0.485,0.25) and 0(2) at (0.224,0.304,0.039), we 
calculated the dipole fields in the unit cell to search for possible muon sites in a locus of 
about 1 A from either the 0(1) or 0(2) sites. Here the atomic coordinates are given in 
fractions of the lattice spacings (5.542 A, 5.732 A, 7.783 A). The calculated dipole fields 
were compared to the measured low-temperature frequencies in LaMnOs, Lao.67Cao.33Mn03 
and CaMnOs. In LaMnOs there are two observed low-temperature frequencies, z/^i = 84.90 
MHz and z/^2 = 128.8 MHz, corresponding to local magnetic fields of 1.6 kOe/fiB and 2.4 
kOe//i£, respectively. We find two muon sites which are consistent with these data. The 
first site [Mu(l), at (0.06,0.71,0.25)] is about 1 A from 0(1) and is quite close to the 
Holzschuh site found for the orthoferrites. The second site [Mu(2), at (0.38,0.34,0.0)] is 
about 1 A from 0(2). The calculated dipole fields in LaMnOs (Mn spin S = 2) for the 
known AFM structure are 1.6 and 2.4 kOe/yUe for Mu(l) and Mu(2), respectively, which 
are in excellent agreement with the measured values. Furthermore, these same two sites 
reproduce the single measured low-temperature frequency in FM La .67Cao.33Mn0 3 , v^ = 75 
MHz, or 1.5 kOe/f^B', i.e., the calculated value of B dip is 1.5 kOe//^B for both of these sites, 
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using S = 2(1 — x) + 1.5x, x = 0.33. Finally, in CaMnOs three muon frequencies are observed: 
V/a. = 14.3 MHz, u^ 2 = 41.2 MHz and v^ = 86.5 MHz, corresponding to 0.35 kOe/// B , 1-01 
kOe/fiB, and 2.13 kOe/fiB, respectively. The calculated dipolar fields yield B dip = 0.27 
kOe/fiB for Mu(2) and B dip = 2.3 kOe/ lib for Mu(l) using S = 1.5, in good agreement with 
Vpx and v^, respectively. Thus, these two muon sites are consistent with all of the measured 
low-temperature muon frequencies in LaMnOs, Lao.67Cao.33Mn03 and CaMnOs except the 
f M 2 frequency in CaMnOs, which corresponds to a metastable position, as discussed below. 
A more definitive identification of the muon sites in these materials would require single 
crystals, in which not only the magnitude but also the direction of the local field could be 
determined. Nevertheless, from the measurements reported here one may strongly infer that 
the positive muon occupies positions within about 1 A of both the 0(1) and 0(2) oxygen 
atoms across the (La,Ca)MnOs series. Because a small amount of local lattice dilation due 
to the electrostatic field of the muon can reduce the calculated dipole fields by about 10%, 
and because we neglect the possible effects of any small transferred hyperfine fields, we 
estimate that our position assignments are uncertain to at least 0.2 A. Although knowledge 
of the muon site is very important for some types of experiments, we point out that the 
interpretation of our primary results, the Mn-ion relaxation rates, does not depend on the 
exact location of the muon in the lattice. 

It is important to know whether the muon is stationary in the lattice, however, because 
rapid muon motion can mimic the relaxation of the Mn spins themselves. The occurence 
of stable muon sites does not apply to CaMnOs, for example. Here three muon frequencies 
were observed, with only the lowest frequency v^ persisting at temperatures as high as T N . 
The amplitude corresponding to v^i also increases monotonically with temperature, as seen 
in the top frame of Fig. ^. This is likely due to a change in the muon site population with 
temperature, as found in the orthoferrites. 

A transfer of polarization between different muon sites occurs as a consequence of the 
thermal activation of the positive muon from a relatively shallow trapping site to a deeper 
one.tJ To affect the observed relative populations the muon hopping time r^ must be 
sufficiently small; e.g., significant population change takes place at a temperature where 
T/jAix^ < 1, with Au; M the frequency difference between the two sitesHI While we do not 
attempt a detailed analysis of the muon motion in CaMnOs, we estimate r^ = 0.003 lis at 
T = 80 K in CaMnOs, assuming Au;^ = 2tt x 50 MHz, corresponding to hopping between 
sites with frequencies u^ and v^i, for example. We note, however, that the 1/Ti linewidth 
is independent of temperature above T/v up to 275 K. This indicates that any postulated 
muon hopping rate out of the most stable site in CaMnOs (corresponding to v^i) is small 
compared to the inverse of the Mn spin correlation time r _1 . An order of magnitude esti- 
mate for r~ x at high temperatures in CaMnOs yields 10 11 — 10 12 s^ 1 . Thus, we deduce that 
3 x 10 8 < r^ 1 < 10 11 - 10 12 s- 1 in CaMnOs for T > T N . 

Having obtained direct evidence for unstable muon sites in CaMnOs, one must consider 
the possibility of muon motion in the other compounds, particularly at the highest measured 
temperatures. Significant muon hopping is ruled out below 170 K in LaMn0 3 , for the 
following reasons. Unlike in the orthoferrites or in CaMnOs, the two signal amplitudes 
which we measure in LaMnOs remain temperature independent below T/v (top frame, Fig. |Tj). 
Thus these data, and the constancy of the spin-lattice relaxation rate above T/v, indicate 
that the muon shows no sign of occupying a metastable site below at least 170 K in LaMnOs. 
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One might inquire whether the falling of the rate X s in Lao.67Cao.33Mn03 above about 267 
K (Fig. 16) could be caused by muon motion. This would require a muon hopping rate 



T, 



1 3> t s 1 , the inverse correlation time of the 'slow' Mn spins. We estimate t s x ph 2 x 10 11 



h ^ 'a 



s 1 , where t s 1 is derived from the measured X s p» 2yus x and the relation A s = 2uj 2 t s , 
using the low-temperature muon frequency u^ = 2tt x 75 MHz. This yields a value for t^ 1 
in Lao.67Cao.33Mn03 which exceeds the upper limit found above for the most stable site in 
CaMn03 at comparable temperatures. (This site, Mu(2), is also found in Lao.67Cao.33Mn03.) 
The disorder induced by doping in (La,Ca)Mn03 should help to stabilize the muon position, 
and it therefore seems unlikely that muon motion is a significant factor in our measurements. 

B. Muon precession frequencies 

The temperature dependence of the v^ reflects the growth of the local magnetization 
below the ordering temperature. We have performed an analysis of the temperature depen- 
dence of the muon frequency below the ordering temperatures in those materials where the 
muon frequency is observable. The temperature dependence of ^(T) can be parameterized 
as follows: 

^(T) = ^ (l-T/T M f (5) 

Fits were carried out in a limited temperature range near the ordering temperature Tm 
(0.80T/U < T < Tm), reflecting the limited range of validity for Eq. (|). (Strictly speaking 
Eq. (|5|) is valid only in the asymptotic critical region of a second order magnetic phase 
transition. Often this expression works well outside of this region, however.) The fits are 
shown as solid lines in the lower frames of Figures [2] and [| for LaMn03 and CaMn03, 
respectively. In LaMn03 the two frequencies were normalized at T = 120 K; one sees 
from the plot in Fig. |2] that the two signals possess the same temperature dependence. We 
obtained (3 = 0.340 ± 0.008 for LaMn0 3 and (3 = 0.53 ± 0.03 for CaMn0 3 ; in the latter case 
only the more stable frequency v^i was included in the fit. The exponent for LaMnOs is 
consistent with measured critical exponents in 3D Heisenberg antiferromagnets,Ej while the 
exponent for CaMnOs yields the mean field value of 1/2. The presence of muon diffusion 
(i.e., change in muon site population) in CaMnOs may render this value unreliable, however. 

The fit for the Lao^Cao.isMnOs material, carried out over the observable temperature 
range (0.8 < T/T c ^S 1.0), is shown in the middle frame of Fig. ||. This yielded T c = 
182.2±1.8 K and /3 = 0.49±0.09. A fit to the La . 67 Ca .33MnO 3 data between T = 220-267 
K yielded T c = 268 ± 0.4 K and (3 = 0.233 ± 0.013.1! 

In each of these cases the fitted frequency curves give a magnetic ordering temperature 
which is consistent with the peak in the spin-lattice relaxation rates arising from critical 
slowing down of the Mn spins. Note that for Lao.67Cao.33Mn03 this is only true for the 
two-exponential fit, not for the stretched exponential fit. This fact is further indication that 
the two-exponential fit is a better representation of the data near Tc for this system. 

It is not clear how to interpret the values of j3 obtained from the FM materials x = 
0.18 and 0.33, for the following reasons. The insulating x = 0.18 material lies close to 
the boundary where canted AFM is observed. Furthermore, the growth of the FM order 
parameter in this compound may be influenced by the change in spin dynamics found below 
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150 K, near where the resistivity shows its anomalous increase. Either of these facts makes 
comparison of /3 with theories of conventional phase transitions dubious. Finally, as discussed 
below, the transition in the x = 0.33 conducting material may proceed via the percolation 
of FM clusters arising from the presence of magnetoelastic polarons, rather than via a more 
conventional FM second-order phase transition. At present there is, therefore, no theory to 
which these exponents can be compared. 



C. Spin- lattice relaxation 

We first discuss the zero- field spin-lattice relaxation rates in LaMnOs and CaMnOs, 
shown in the top frames of Figs. |2] and f|, respectively. In both LaMnOs and CaMnOs G r i x 
is well described by a single exponential A r \ x exp(—t/Ti). The peak in the spin-lattice relax- 
ation rate 1/Ti reflects critical slowing down near T/v, and in each material the temperature 
region of the critical dynamics is relatively narrow. Furthermore, the zero-field transition 
in both materials is quite sharp, as seen by the rapid decrease in the relaxing amplitude 
A r i x below T/v, shown in the middle frames of Figs. ^| and |j. Note that even the relatively 
modest (hbH <C /c_bTv) field H = 3 kOe applied parallel to the muon spin strongly affects 
the critical dynamics. 

Both LaMn03 and CaMn03 display temperature independent spin-lattice relaxation over 
an extended temperature range above T/v- This is characteristic of systems fluctuating at 
their exchange frequency u e . In this instance one hast3 



1/Ti = V2ttS(S + l)A7(3/rw e ), (6) 

where S is the Mn atomic spin (S = 2 for LaMnOs and S = 3/2 for CaMnOs) and A is the 
hyperfine coupling constant between the muon spin and the spin S. The exchange frequency 



is related to the exchange integral J via the expression hu e = Jj2zS(S + l)/3, where z is 
the number of nearest neighbor Mn spins to a given Mn spin. 

We now use Eq. (|6]) to estimate the exchange frequency for LaMnOs and CaMnOs. To 
do so one may obtain A from the value of the low temperature muon frequency (27rz/ M = u^) 
or frequencies. If there is a single frequency then A = hu^/S. When there is more than one 
muon frequency the relaxation function takes the form G rlx = £j(v4 r i x )jexp(— t/T u ), where 
the index i runs from 1 to 3 in CaMnOs, for example. In principle, this implies that one 
should observe a multi-component or stretched-exponential relaxation function, which is in 
contradiction to the single-exponential function observed in both LaMnOs and CaMnOs. 
However, if the individual relaxation times Tu are all long compared to the characteristic 
time of the experiment (of the order of the muon lifetime), then 

G r i x = £j(A rlx )jexp(-t/Tij) w v4 r i x exp(-£/Ti), (7) 

where 1/Ti = E^Arix^T-jJ /.A r ix- in this limit a single exponential relaxation function is 
recovered, with a relaxation rate which is an appropriately weighted sum of rates 1/Tij. The 
observed relaxation times T\ in both LaMnOs and CaMnOs are relatively long (« 10 — 20/xs), 
and thus the individual relaxation times Tu are also presumably long. This then explains 
the observation of an exponential relaxation function in each case. 
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In the motional narrowing limit ,113 appropriate for relaxation above T/v, 1/Tu = 2oj 2 i Ti, 
so that A 2 m Yi^A^^Tiuj^/ S) 2 / A T \ X . Here r is the Mn-ion correlation time. In LaMn03 
the observed low temperature frequencies are uo^i = 2tt x 84.9 MHz and uo^ = 2tt x 128.8 
MHz. These values yield an exchange frequency u e = 6.33 THz. corresponding to an ex- 
change integral J = 0.85 meV. Using the mean field expression^ T/v = JzS(S + l)/3 one 
obtains T/v = 119 K, in good agreement with the measured Neel temperature of 139 K. The 
situation is not so sanguine in CaMnC>3, however. Here the low temperature frequencies are 
cu^a = 2tt x 14.3 MHz, u^ = 2vr x 41.2 MHz and uj^ = 2vr x 86.5 MHz, yielding uj e = 2.8 THz, 
J = 0.47 meV and T/v = 42 K, lower by a factor of about 3 than the measured T/v- However, 
if one uses only the lowest frequency v^i for CaMn0 3 (corresponding to the 'stable' muon 
site), one obtains uo e = 0.10 THz, J = 0.017 meV and T/v = 1.5K, an unacceptably low value 
even given the crudeness of the approximations used. This seems to be an indication that 
in CaMnC>3 muon motion makes determination of the effective coupling constant A ambigu- 
ous. Alternatively, the mean field approximations used to estimate T/v may not be valid in 
CaMnC>3, despite the finding of a mean- field exponent for the temperature dependence of 
the lowest /xSR frequency Vp\. 

As the amount of Ca doping increases from zero in La 1 _ rr Ca z Mn03 the width of the 
peak in 1/Ti near the magnetic transition grows larger. This is seen already for x = 0.06 
(Fig. ^j, where the width in 1/Ti extends over a considerably larger temperature region 
than does the fall in the relaxing amplitude A-ix- Recall that a reduction in A r \ x signals the 
onset of magnetic order, even though no muon frequencies were observed in this material. 
The changing value of the exponent K also extends over a wide temperature regime around 
T/v, indicating that the spin dynamics are becoming inhomogeneous near T/v. This trend, 
a general decline in the value of K at the magnetic ordering temperature, is seen in all 
of the Ca-doped samples except the end compound CaMnC>3, in which simple exponential 
relaxation is observed at all temperatures. Furthermore, the 1/Ti relaxation rates in the 
x = 0.06 and x = 0.33 materials tend to remain relatively large at temperatures significantly 
below their magnetic ordering temperatures, indicating that there is an anomalously slow 
(large r) relaxation process which persists in the ordered state. Values of K <C 1, combined 
with the real-space nature of the muon probe, suggest that the relaxation rates are therefore 
spatially distributed, as discussed in more detail below. 

The FM x = 0.18 compound and the AFM x = 0.67 compound each exhibit an uncharac- 
teristic increase in their spin-lattice relaxation rates well below their ordering temperatures. 
This change is marked by a broad maximum near T « 110 K for x — 0.18 (Tq = 182 K) and 
T « 85 K for x = 0.67 (T/v — 150 K). It is tempting to attribute the loss of the oscillating 
component and the gentle rise in the relaxation rates near 110 K in Lao^Cao.isMnOs to 
the influence of 'charge ordering' (CO) reported by otherso at T rj 70 K. It is known, for 
example, that the dynamical spin correlations in (Bi,Ca)Mn03 change from ferromagnetic 
above the CO temperature T co to antiferromagnetic below T co .E3 We note, however, that 
CO observed in La .33Cao.67Mn0 3 occurs at temperatures above T N , as opposed to below 
Tc in Lao.82Cao.i8Mn03. Because /xSR cannot determine the type of magnetic correlations 
present in these materials (in general, this requires a momentum-dependent probe), we can- 
not definitively explain the change in spin dynamics below the ordering temperatures of 
the x = 0.18 and x = 0.67 materials. We do note the possible correlation with CO in 
both, however. This could be the subject of further investigation by neutron scattering, for 
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example. 

The two FM materials under study, Lao^Cao.isMnOs and Lao.67Ca .33Mn0 3 , exhibit a 
bi-modal distribution of fluctuation rates in perhaps the most interesting part of the phase 
diagram. The two-exponential fits to G v \ x [Eq. (f|)] yielded relaxation rates and relative 
amplitudes A/, X s , Af and A s . We have labeled these quantities 'f and 's' to indicate 
'fast' and 'slow' Mn-ion fluctuation rates r _1 because, in general, the local muon relaxation 
rate A ex 7^Z) g \^B(q)\ 2 r(q), as mentioned in Section IA. Here \SB(q)\ is the amplitude of 
the local fluctuating field and r(q) is the Mn-ion correlation time. The relaxation rates Xf 
and A s display very different magnitudes (Figs. [15] - [161): in each material X s ~ 10 — 40 x A/. 



Using the relation A = luipj, one may derive Mn correlation times for Lao.67Cao.33Mn03 
at T = 250 K, for example. One obtains r/ ~ 10~ 13 s and t s = 5 x 10~ 12 s, assuming 
that uj^ at the lowest temperatures is representative of the /i + -Mn coupling for both the 
slow and fast components. It is noteworthy that only Xf shows evidence of 'critical slowing 
down' at Tc. In the FM conductor Lao.67Cao.33Mn03, Af increases below Tc (and A s 
decreases), while in the insulating material Lao^Cao.isMnOs the fast and slow amplitudes 



are relatively temperature independent below Tc (Fig. |17[). At Tc, Af is also a somewhat 
larger in Lao.82Ca .i8Mn0 3 than in La .67Cao.33Mn0 3 . 

This analysis is consistent with two different spatially distinguishable regions in the sam- 
ple, characterized in our measurements by very different relaxation rates and temperature- 
dependent volumes. More specifically, although the Xf component may exist throughout 
the sample volume, the X s component can only exist in spatially separated regions; other- 
wise only the slow Mn fluctuation rate (which creates the largest /iSR rate and corresponds 
to X s ) would be observed. The muon samples these separate regions locally through its 
relatively short-ranged dipolar coupling to the Mn spins (matrix element ~ r~ 6 ). It is im- 
portant to note that neither /iSR relaxation component corresponds to pure LaMnOs or 
CaMn03, whose relaxation rates are independent of temperature above 140 K with magni- 
tudes ~ 0.10 — 0.05 fis^ 1 , as discussed earlier. This argues against large-scale separation of 
Mn 3+ and Mn 4+ ions. Furthermore, we note that because of the large spin wave stiffness 
constant observed^ near Tc in these FM materials the muon spin is not significantly relaxed 
by spin waves, as pointed outt3 previously. Thus, these FM materials possess at least three 
types of fluctuations, the FM spin waves observed in neutron scattering, and the excitations 
associated with the Xf and X s relaxation rates reported here. 

We postulate that the two relaxation components observed here constitute the spin signa- 
ture of magnetoelastic polarons as the system becomes FM. As discussed in the introduction 
to this paper, these polarons are thought to be formed from FM polarized Mn spins which 
surround a region of local lattice distortion, and presumably originate from the combined 
DE and JT interactions. These polaronic effects, together with the random replacement of 
La atoms by Ca, produce considerable disorder. Our measured A/ component shows a peak 
in the relaxation rate at Tc, and we associate this relaxation with overdamped spin waves, 
characteristic of a disordered FM. That is, the damping of the spin waves in this scenario 
would be caused by a loss of translational symmetry arising from inherent disorder. 

As mentioned in the introduction, a broad peak in the neutron scattering intensity cen- 
tered around zero energy transfer has been observed in FM Lao.67Cao.33MnO3.li3 This peak 
coexists with Mn spin waves near T c , and persists to temperatures as low as T/T c = 0.7, 
consistent with the finite /zSR rate we measure in Lao.67Cao.33Mn03 below Tc- The neutron 
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linewidth T obeys a diffusive relaxation law, T = Dq 2 , where D rj 30 mev A 2 . The length 
scale for this diffusion was found to be 12 A from small angle scattering.!^ We can compare 
this linewidth T with A/ using the relatione^ 

X f = 2u 2 f(S(q, t) ■ S(-q, 0)) exp(iwt)dt. (8) 

Here (S(q,t) ■ S(— q, 0)) is the transverse spin-spin correlation function given by 

(S(g, t) ■ S(-g, 0)> = k B T X (q) exp(-Dq 2 t)/(gfi b ) 2 , (9) 

where x(l) is the static g-dependent susceptibility k B is Boltzman's constant and \xb is the 
Bohr magneton. We take an Ornstein-Zernike form for x{o)i given byc3 

X (q) = (g^ b ) 2 S(S + l)/(3k B T(q 2 + T 2 )r 2 ), (10) 

where r\ = a 2 /6 for a cubic lattice with lattice constant a, and £ is the coherence length 
for the diffusive spin correlations. Evaluating these expressions for £ = 12 A yields Xf = 
0.18 /xs -1 , in excellent agreement with our /iSR measurements of A/ at T = Tq (see Fig. 
16|) . Thus, we associate our observed Aj with the diffusive relaxation mode observed in 
quasi-elastic neutron scattering, and attribute the relaxation to overdamped spin waves. 

The existence of overdamped FM spin waves is naturally associated with a relatively 
high level of charge mobility because the FM in these CMR compounds is induced by the 
hopping of holes via the DE interaction. By contrast, the X s component corresponds to a 
rather long correlation time of ~ 10~ n s. This order-of-magnitude is probably inconsistent 
with electronic relaxation of single Mn 3+ or Mn 4+ spins. Other possible mechanisms might 
include phonon induced spin relaxation (via magnetoelastic coupling) or the slow overturning 
of spins associated with local polaron hopping in regions of the sample where spin and charge 
motion are frustrated by more extreme local lattice distortions, that is, regions which are 
relatively insulating. 

The possiblity that the X s component is associated with relatively insulating regions 
of the sample, while the A/ component is associated with relatively conducting regions, is 
qualitatively consistent with the temperature dependence of the volume fractions A s and 
Af in the two FM materials. In the FM conductor Lao.67Cao.33Mn03 the magnitude of 
Af increases (and A s decreases) as the temperature is lowered below Tq and the material 
becomes more conducting. By contrast, in the insulating FM Lao^Cao.isMnOa the volume 
fractions Af and A s are independent of temperature. This suggests that even significantly 
below Tq in conducting Lao.67Cao.33Mn03 there exist relatively non-conducting regions which 
are the remnants of the inhomogeneity produced by isolated polarons existing above Tc in 
these materials. These less conducting regions could be caused by an uneven distribution 
of Ca atoms, or by more extreme local lattice distortions which inhibit charge motion. As 
the temperature is lowered below Tc in Lao.6?Cao.33Mn03 these less conducting regions are 
gradually absorbed into the larger conducting matrix of the material. This apparently does 
not occur in the insulating FM Lao^Cao.isMnOs, however, where the volume fractions are 
frozen below Tc- 

Such a postulated scenario is consistent with our measurements, and with charge trans- 
port and structural measurements, as follows. In the two-fluid model of the resistivitynJ 
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growth of the conducting charge fraction in Lao.67Cao.33Mn03 occurs rapidly as the mag- 
netization grows, reaching essentially 100 % just below Tq. This does not contradict the 
above picture from the /iSR data (in which slowly fluctuating spins in relatively less metal- 
lic regions persist far below T c ) because when a conducting path is reached (at To) the 
resistivity is shorted out, even though a considerable volume fraction may still be relatively 
non-conducting. Thus the two-fluid model of the resistivity is similar to the picture pre- 
sented here, with this one exception. The persistence of an inhomogeneous ground state 
below Tc is reinforced by local probes of the lattice structure,t!rE3 which are consistent 
with the gradual loss of structural inhomogeneity below Tq, similar to that reflected in the 
spin-lattice relaxation rates reported here in the conducting FM Lao.67Cao.33Mn0 3 . 



IV. SUMMARY 

In this paper we present /iSR data for a series of (La,Ca)Mn0 3 compounds, span- 
ning the range from the AFMs LaMn03 and CaMn03, to the FMs Lao^Cao.isMnOs and 
La .67Cao.33Mn03, to the charge-ordered material Lao.33Cao.67Mn03. The /iSR technique 
allows one to obtain the temperature dependence of both the local static magnetic order pa- 
rameters and the Mn spin-lattice-relaxation rates r _1 , where generally 10~ 4 < r _1 < 10~ 13 
s _1 . These measurements, therefore, provide a necessary complement to other spin probes, 
such as neutron scattering and nuclear magnetic resonance. 

From the magnitude of the measured local magnetic fields we have identified two plausible 
stable muon sites near the 0(1) and 0(2) atoms in these materials. Except for CaMn03, 
where evidence for metastable muon sites is found, we conclude that rapid muon diffusion 
is not a factor in our measurements of the Mn spin dynamics. We find that the spin 
dynamics near the AFM and FM magnetic ordering temperatures of these compounds is 
quite sensitive to small applied magnetic fields, however. This may be a consequence of the 
FM double-exchange interaction, which is enhanced by the application of a uniform field. 
From the magnitude of the zero-field ^SR rate above Tjy in the AFM insulating parent 
compound LaMn03 we derive an exchange integral J = 0.83 meV, which is consistent 
with the mean field expression for Tjy. Our principal finding is that the addition of Ca 
to LaMn0 3 (or La to CaMn0 3 ) results in an inhomogeneous spectrum of zero- field /iSR 
rates below the magnetic ordering temperatures in these materials. In the FM compounds 
Lao.82Cao.i8Mn03 and Lao.67Cao.33Mn03 we are able to resolve this inhomogeneity into a 
bi-modal distribution of 'fast' and 'slow' /zSR rates, corresponding to Mn spin fluctuation 
rates which differ by about a factor of at least 30. Because /zSR is a local probe this implies 
that the 'slow' Mn fluctuations must exist in distinct spatial regions of the sample. A 
physical picture is hypothesized for these FM phases in which the fast Mn rates are due 
to overdamped spin waves characteristic of a disordered FM, and the slower Mn relaxation 
rates derive from distinct, relatively insulating regions in the sample. These results present 
direct evidence for microscopic inhomogeneity in the spin channel of the FM manganites, 
suggesting that a proper description of these materials probably includes more than just the 
simple JT and DE interactions. 
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FIG. 1. Temperature dependence of the fiSR oscillating amplitude (top), precession frequencies 
(middle) and fractional linewidth (bottom) in LaMn03 from fitting to Eq. (g). 

FIG. 2. Temperature dependence of the /uSR spin-lattice relaxation rate (top), relaxing am- 
plitude (middle) and muon precession frequencies v^x and v^ (bottom)in LaMn03. v^i has been 
normalized to Vn\ at 120 K. The solid curve in the bottom frame is from a fit of the precession 
frequencies to Eq. (|5|). 

FIG. 3. Temperature dependence of the //SR oscillating amplitude (top), precession frequencies 
(middle) and fractional linewidth (bottom) in CaMn03 from fitting to Eq. (g). 

FIG. 4. Temperature dependence of the /iSR spin-lattice relaxation rate (top), relaxing ampli- 
tude (middle) and muon precession frequency u^i (bottom) in CaMn03. The solid curve in the 
bottom frame is from a fit of v^i to Eq. (||). 

FIG. 5. Temperature dependence of the resistivity for Lai^^Ca^MnOs, x = 0.0,0.06,0.18 and 
0.33. Comparable data for the x = 0.67 sample are given in S-W. Cheong and C. H. Chen, Colossal 
Magnetoresistance and Related Properties, ed. B. Raveau and C. N. R. Rao (World Scientific), to 
be published. 

FIG. 6. Temperature dependence of the ^SR relaxing amplitude (top), exponent (middle) and 
spin-lattice relaxation rate (bottom) in Lao.94Cao.o6Mn03 obtained from a fit to Eq. (y). The filled 
symbols are for zero applied field, and the open circles for an applied field of 2.5 kOe. 

FIG. 7. Temperature dependence of the /iSR relaxing amplitude (top), exponent (middle) and 
spin-lattice relaxation rate (bottom)in Lao.82Cao.isMn03 obtained from a fit to Eq. (||). 

FIG. 8. Temperature dependence of the /uSR fractional linewidth (top), precession frequency 
(middle) and spin-lattice relaxation rate (bottom)in Lao.82Cao.isMn03. The data in the top and 
middle frames are from a fit to Eq. (||). The solid line in the middle frame is from a fit of v^ to 
Eq. (§). 

FIG. 9. Temperature dependence of the ^SR fractional linewidth (top), precession frequency 
(middle) and spin-lattice relaxation rate (bottom) in Lao.67Cao.33Mn03 obtained from a fit to 
Eq. (^). The inset in the middle frame shows a fit of v^ to Eq. (|B|). 

FIG. 10. Temperature dependence of the /iSR relaxing amplitude (top), spin-lattice relaxation 
rate (middle) and exponent (bottom) for zero applied field (filled symbols) in Lao.67Cao.33Mn03 
obtained from fits to Eq. (pi). The open symbols in the middle frame show data in applied magnetic 
fields of 1 kOe and 3 kOe. 

FIG. 11. Temperature dependence of the /iSR relaxing amplitude (top), spin-lattice relaxation 
rate (middle) and exponent (bottom) for zero applied field in Lao.33Cao.67Mn03 obtained from fits 
to Eq. (§. 

FIG. 12. (a) /iSR relaxation function G r i x (t) at T = 270 K. The curves show best fits using 
stretched exponential from Eq. (||) [1/Ti = 0.219(7) /is -1 , K = 0.63(3)] and exponential (K = 1) 
functions, (b) A two-exponential least squares fit using Eq. @, with Af = 0.60(5), A s = 0.40(5), 
X f = 0.106(5) [is" 1 and A s = 0.93(11) fis' 1 . 
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FIG. 13. The measured ^SR relaxation function G r \ x (t) at T = 270 K (open squares), normal- 
ized such that G r i x (0) = 1. The lines show calculated stretched-exponential (lower dotted line) 
and two-exponential (upper dashed line) functions using the fit parameters obtained from the fits 
shown in Figure |i"2]. A stretched exponential A r \ x exp(— t/T\) K is linear with slope (1 — K) on a 
plot of this type. 



FIG. 14. Time dependence of the ^SR relaxing asymmetry G r \ x (open squares) for zero applied 
field in Lao.67Cao.33Mn03 showing the curves Xf (slowly relaxing curve corresponding to faster Mn 
spins) and X s (rapidly relaxing curve corresponding to more slowly relaxing Mn spins) for a two 
exponential fit. The curve for X s has been offset from the data for clarity. 



FIG. 15. Temperature dependence of the ^SR fast (Xf, top) and slow (A s , bottom) spin-lattice 
relaxation rates for zero applied field in Lao.82Cao.i8Mn03 obtained from a fit to Eq. (y) [closed 
symbols, exponent K ~ 1 in Eq. (y)] and Eq. (||) [open symbols]. 



FIG. 16. Temperature dependence of the //SR fast (Xf, top) and slow (A s , bottom) spin-lattice 
relaxation rates for zero applied field in Lao.67Cao.33Mn03 obtained from a fit to Eq. (^) [closed 
symbols, exponent K ~ 1 in Eq. (||)] and Eq. (||) [open symbols]. 



FIG. 17. Temperature dependence of the /xSR fast (Af, top) and slow (A s , bottom) amplitudes 
for zero applied field in Lao.67Cao.33Mn03 and Lao.82Cao.i8Mn03 obtained from a fit to Eq. (|3|) 
[closed symbols, exponent K ~ 1 in Eq. @] and Eq. (||) [open symbols]. 
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